Transverse Expansion and High p T Azimuthal Asymmetry at RHIC 
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Rapid 3+1D transverse plus Bjorken collective expansion in A + A collisions at ultra-relativistic 
energies is shown to reduce substantially the azimuthal asymmetry resulting from jet quenching. 
While the azimuthal asymmetry in non-central collisions, V2(pr > 2 GeV/c) ~ 0.15 reported by 
STAR at RHIC, can be accounted for by spatially anisotropic jet energy loss through a 1+1D 
expanding gluon plasma with dN 3 jdy ~ 1000, we show that if rapid transverse collective expansion 
of the plasma is assumed, then the asymmetry due to jet quenching may be reduced below the 
observed level. Possible implications of this effect are discussed. 

PACS numbers: 12.38. Mh; 24.85. +p; 25.75.-q 



I. INTRODUCTION 



Preliminary STAR data from RHIC §-§ suggest a 
large transverse azimuthal asymmetry i>2 ~ 0.10 — 0.15 
at moderate high 2 < pr < 5 GeV/c. In Rcfs. [|]|]] 
some of us proposed that the saturation of V2(pt) could 
be interpreted as a manifestation of asymmetric jet en- 
ergy loss that arises due to the different jet propaga- 
tion path lengths in non-central (b > fm) collisions 
of heavy ions. We applied the GLV [7|-|ir^] energy loss 
formalism to compute the asymmetry and found that un- 
der the approximation that the initial spatial anisotropy 
of the plasma (generated by intersecting sharp cylin- 
ders) did not change with time the observed high px «2 
could be explained by asymmetric jet quenching. The 
initial gluon rapidity density when extrapolated to cen- 
tral collisions was estimated to be dN g /dy(h — fm) ~ 
1000. The implied large gluon density, Pg(r, b) » 
dN 9 /dy(b)/(rAj_(b)), was also found to be consistent 
with the azimuthally averaged jet quenchi ng p attern of 
7T° at px > 2 GeV/c reported by PHENIX tfutfL^. 

In Ref. D the non-perturbative low pt component of 
the hadronic spectra was parametrized in terms of hy- 
drodynamic solutions QJl^l that fit well the observed 
Pt < 2 GeV/c flavor dependent inclusive distributions 
for pions and nucleons We ignored in Ref. ||, 

however, the additional rapid dilution and time depen- 
dent spatial asymmetry of the plasma that is predicted by 
hydrodynamic transverse flow. Boosted thermal model 
fits |i~7|-pl|] to the low data suggest that the collective 
transverse flow velocities may be quite high at RHIC with 
vt ~ 0.6. The purpose of this letter is to compute the 
effect of such 3+1D expansion on the azimuthal asym- 
metry of the jet quenching pattern. We find that while 
in extreme scenarios, assuming instantaneous transverse 
expansion of the medium, the geometric anisotropy is 
strongly reduced, the overall magnitude of the mean en- 
ergy loss AE is much less sensitive for flow velocities 
up to vt — 0.8. This can be verified for a variety of 
transverse density profiles as well as a full hydrodynamic 
calculation 



22 23f| . We also speculate about differ- 



baryon production that could lead to different flavor de- 
pendent high pt anisotropy behavior of hadron yields. 



II. ENERGY LOSS OF JETS IN TRANSVERSELY 
EXPANDING MEDIUM 

The explicit solution to the problem of energy loss of a 
hard ultra-relativistic jet (c = 1, z = r) produced inside 
a hot and dense medium at position xo and propagating 
in direction v was obtained in jq| to all orders in opacity 
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dr ct(t)p(xo + v(t - To), r) 



(1) 



It was shown that for gluon plasma of transverse thick- 
ness on the order of the nuclear size and density p < 
50 fm -3 (relevant at the currently available heavy ion 
collider energies) the dominant contribution comes from 
the first order in the opacity expansion |Q,||j24|. In the 
case of 1+1D Bjorken longitudinal expansion with ini- 
tial plasma density po = p(zq = r ) and formation time 
z = t , i.e. 
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it is possible to obtain a closed form analytic formula || 
(under the strong asymptotic no kinematic bounds as- 
sumption) for the energy loss due to the dominant first 
order term For a hard jet penetrating the quark- 

gluon plasma 
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where x ~ uj/E is the momentum fraction of the radi- 
ated gluon and the formation physics function f(Z(x, t)) 
is defined in [El to be 
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du 
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ent non-perturbative effects on moderate pt meson and 



With Z(x,t) = (t - t )(i 2 (t)/2xE being the local 
formation physics parameter, two simple analytic lim- 
its of Eq. (Q) can be obtained. For x x c = 
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p(t ) 2 t 3 L 1 -^/2E = Lfi 2 (L)/2E, in which case the 
formation length is large compared to the size of the 
medium, the small Z(x, r) limit applies leading to f(Z) w 
irZ/2. The interference pattern along the gluon path be- 
comes important and accounts for the the non-trivial de- 
pendence of the energy loss on L. When x <C x c , i.e. the 
formation length is small compared to the plasma thick- 
ness, one gets f(Z) w I02Z. In the x ^> x c limit ^ the 
radiative spectrum Eq. (J3|) becomes 
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dx 2(2 -a) A(t ) x 

The mean energy loss (to first order in x) integrates to 
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The logarithmic enhancement with energy comes from 
the x c < x < 1 region ffijM. In the case of sufficiently 
large jet energies (E — > oo) this term dominates. For 
parton energies < 20 GeV, however, corrections to this 
leading log l/x c expression that can be exactly evaluated 
numerically from the GLV expression and are found to be 
comparable in size. The effective AE/E in this energy 
range is approximately constant pj| . It was also recently 
shown that the effect of gluon absorption is negligible for 
jet energies E > 5p ~ 3 GeV @. 

In the above expression, transverse expansion can only 
be very crudely modeled taking a > 1. To derive an im- 
proved analytic expression taking transverse flow into ac- 
count, we consider next an asymmetric expanding sharp 
elliptic density profile the surface of which is defined by 
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The area of this elliptic transverse profile increases with 
time , r, as 



A±(t) = n(R x + v x T)(R y + v v t) 



(8) 



Consider the plasma density seen by a jet in direction 
(j)o starting from xo = (xo,yo) inside the ellipse with a 
specified initial density ropo = ^/(^RxRy) dN 9 /dy. The 
density along its path is 
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Analytic expressions can be obtained only for asymp- 
totic jet energies when the kinematic boundaries can be 
ignored. In that case Eq. (||) reduces to 
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where we kept only the dominant log l/x c , Z < 1 regime. 
By integrating over x the total energy loss is for this den- 
sity profile azimuthally dependent: 
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which is a linear (mod log) weighed line integral over the 
asymmetric density. This integral is still convergent in 
spite of the singularity of the density at r = 0. 

We consider in more detail the analytically tractable 
case of a sharp expanding elliptic cylinder. We 
approximate the assumed 4>o independent screening 
p(r) w gT(r) = 2(p(r)/2) 1 / 3 since g ~ 2 and p = 
(16C(3)/tt 2 )T 3 ~ 2T 3 for gluon plasma. We define t(4> q ) 
as the escape time to reach the expanding elliptic surface 
from an initial point xq in the azimuthal direction </>q: 
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We take to(4>o) = 2r(</> ) (/°( T (0o))/2) 2 ^ 3 to estimate an 
upper bound on the logarithmic enhancement factor. 
Performing the remaining integrals one gets: 

a nfi)/, \ 9C R a 3 s dN 9 E 
4 dy uj(<Po) 
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where a x — v x /R x ,a v = v y /R y . This expression is the 
central result of this paper and provides a simple analytic 
generalization that interpolates between pure Bjorken 
1+1D expansion for small a XtV r, and 3+1D expansion 
at large a XiV T. 

In the special case of pure Bjorken (longitudinal) ex- 
pansion with v x — v y = 
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In this case, the energy loss depends linearly on r((f>). 

Another special case is azimuthally isotropic expansion 
with R x = R y = R and v x = v y = vt- Taking also the 
longitudinal Bjorken expansion into account leads in this 
case to 
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We note that for a jet originating near the center of the 
medium and fully penetrating the plasma the enhanced 
escape time due to expansion r = R/(l — vt) compen- 
sates for the 1/(1 4- vtt(4>o)/R) dilution factor. There- 
fore, in this isotropic case, the extra dilution due to trans- 
verse expansion has in fact no effect of the total energy 
loss 



AE[^(b = fin) w A£^(b = fm) 
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modulo logarithmic factors. We note that in performing 
the line integral Eq. (|l^), the logarithm dependence on 
the cut-off, uj(4>o) was neglected. For our numerical es- 
timates we approximate w(0o) by its average. Our first 
important conclusion is that the inclusive azimuthally 
averaged jet quenching pattern in central collisions is ap- 
proximately independent of transverse expansion. 

It is important to verify Eq. ( |l6| ) for more realistic den- 
sity profiles since it implies that the overall suppression of 
the high pr particle spectra is not affected by transverse 
flow. We have checked numerically that this approximate 
line integral Eq. ( [TT| ) independence of the transverse ex- 
pansion velocity vt holds for transverse density profiles 
with diffuse edges, namely a Gaussian profile 



p(r, (j} ;x ,y ) 



1 



exp 



~ ~dy~ [t) v{R e J f + v xT )(R e y ff + VyT) 
(x + tcos0 o ) 2 (t/o + Tsin0 o ) 2 



(R e J f + v xT f (R e y ff + VyTf 



(17) 



and a curious (0o asymmetric) exponential profile 
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One of the most interesting cases arises as a solution 
of Vlasov equation for a free-streaming massive partons. 
We consider the ideal Bjorken case with perfect corre- 
lation between kinetic and space-time rapidities: y = r) 
with y = tanh -1 (p z /po), V = tanh _1 (z/t) and trans- 
verse propagation at mid-rapidity (r = t). Given initial 
transverse density po(x) and with azimuthally isotropic 
angular distribution of plasma particles with a fixed mag- 
nitude of velocity, at time r one gets 
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For Gaussian initial density the evolution is actually 
given by 
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This solution is rather different from Eq. (|lj]) because 
the density expands as a thick shell rather than homo- 
geneously. The same observation applies for the sharp 
elliptic profile Eq. (^]). The general line integral that 
is needed assuming free Vlasov expansion but averaging 
over initial coordinates and the angle of the internal par- 
ticle motion for fixed azimuthally isotropic vt is 

L(n) = / cPx T po(xo,T ) / dr rp(5c Q + nr,r) 

J J m 



t / d 2 x T po(xo,T ) 



dft r 



2tt 



dr po(x + (n - vt)t, t ) . (21) 



The explicit averaging over the jet production point x i n 
this equation can only be performed numerically. How- 
ever for a jet originating near the center of the medium 
one can derive expression similar to Eq. (|l5[) 
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The equivalent escape factor (enhanced propagation 
length modulated by the dilution of the medium) is given 
in this Vlasov scenario by 
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where K(x) is the complete elliptic integral of first kind. 
The escape factor stays within ±10% of R e ^ for flow 
velocities up to vt — 0.8. As vt — > 1, T(0\))vias. diverges 
only logarithmically. In case of non-central collisions, 
b / 0, and averaging over the jet production points, 
(A^W) . , the validity of Eqs. can be confirmed 

numerically. 

In non-central collisions, the azimuthal asymmetry of 
the energy loss can be expanded in Fourier series and 
characterized as 



A£^(0) = AE(l + 28 2 (E) cos20+---) 
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From Eq. (filf), we compute numerically the azimuthally 
averaged mean AE and the energy loss asymmetry S 2 (E) 
averaging over the uniform elliptic density, Eq. (|7|), of 
the jet formation points. We consider only the specific 
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case of a E = 10 GeV gluon jet to illustrates the de- 
pendence of these quantities on the mean expansion ve- 
locity. The spatial anisotropy is fixed by R x = 2.5 fm 
and R y — 4.7 fm that approximately corresponds to 
the same second moment of the spatial anisotropy of 
the sharp intersecting cylinders considered in Ref. |j| 
for b = 7 fm. The gluon plasma rapidity density that 
fixes the magnitude of the energy loss was assumed to be 
dN g /dy (b = 7 fm) = 325 with a s = 0.2 for illustration. 

We calculate the energy loss and its asymmetry also 
by using a full hydrodynamic calculation of Ref. [flii . 
In this case we use the parametrization eBC of Ref |1J] 
to initialize the system and treat gluon number as con- 
served current to calculate the density evolution needed 
in the line integral Eq. ([H]). We again average over the 
jet formation points but in this case their density is not 
constant but given by the number of binary collisions per 
unit area as in the Woods-Saxon geometry of Ref || . 

We plot in Figs. 1 and 2 the dependence of AE and 
62 (E = 10 GeV) on the mean transverse expansion veloc- 
ities, fixing the velocity asymmetry from hydrodynamic 
results to be v x = + v a ) and v y — Ut(1 — v a ). 

Boosted thermal sources simulations and hydrodynamic 
calculations indicate that vt ~ 0.5 — 0.6 (|Jl^,|2^] and 
v a 0.05 — 0.10. The analytic computation and the 
hydro simulation in Fig. 1 were normalized to the same 
mean energy loss for the vt — point. We see that AE 
is weakly dependent on transverse flow for the schematic 
but analytically tractable uniform expanding elliptic den- 
sity and we have obtained similar results for a wide vari- 
ety of transverse profiles Eqs. ( ^7|Jl^j20| ) . Full hydrody- 
namic calculation of the cf> and xrj averaged line integral 
Eq. ([ll]) is also consistent with approximately constant 
mean energy loss within ±5% for both central and semi- 
central collisions. 

However, the azimuthal asymmetry of the energy loss 
is strongly reduced for realistic hydrodynamic flow ve- 
locities. This implies a much smaller vi at high pt than 
obtained in Ref. Q where transverse expansion was not 
considered. Hydrodynamic calculations with transverse 
expansion predict similar reduction of the averaged line 
integral Eq. (pi]). We note that factor of 1.5 - 2 differ- 
ence in 62(E) without transverse flow in Fig. 2 reflects 
the difference between the sharp cylinder geometry and 
the Wood-Saxon geometry as shown in Ref. M. 
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FIG. 1. The azimuthally averaged energy loss for a 10 GeV 
gluon jet propagating through a 3+1D elliptic expanding 
plasma is plotted as a function of the mean transverse flow 
velocity, vt- The initial profile is a homogeneous ellipse ap- 
proximating Au + Au at b = 7 fm impact parameter. Trans- 
verse velocity asymmetry is fixed to be u a — 0.1. 1+1D and 
3+1D hydro calculations are shown. 
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FIG. 2. The second harmonic of the jet energy loss for a 
10 GeV gluon propagating through a 3+ ID elliptic expand- 
ing plasma as a function of the mean transverse flow veloc- 
ity, vt, is shown. The initial profile is chosen as in Fig. 1. 
We study 3 different transverse velocity asymmetries, i.e. 
u a = 0.0, 0.1, 0.2. Hydrodynamic 1+1D and 3+1D calcula- 
tions using diffuse Wood-Saxon geometry are also included. 



The error bars in the hydro calculation are our esti- 
mate for uncertainties in terminating the line integral 
Eq. ([ll]) at freeze-out and the behavior of density during 
the mixed phase and hadronization. The correspondence 
of the expansion velocity vt in hydro and schematic 
model is also ambiguous since in hydro vt is the aver- 
age velocity on the freeze-out surface, not the velocity at 
the edge of the system as in Eq. (ph. 



III. HIGH P T AZIMUTHAL ASYMMETRY 
SCENARIOS 

The results of the previous section suggest that if 
strong transverse expansion of the plasma can be con- 
firmed (and other scenarios including soft string fragmen- 
tation, baryon junctions |l9|;^0|,^7| and classical Yang- 
Mills pa] evolution rejected) then the high transverse 
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momentum asymmetry predicted by energy loss calcu- 
lations may be reduced (by a factor of 2 to 4). This 
would leave open the understanding of the present STAR 
data in terms of the eikonal picture ||. Interesting 
new possibilities arise to explain the moderate to high 
Pt dependence of the transverse azimuthal asymmetry 
V2(pr > 2 GeV/c). We suggest that the key to con- 
verging to the correct dynamical picture is the flavor de- 
pendence of both the asymmetry and the px differential 
particle multiplicities. We discuss three possible scenar- 
ios for high pt hadron production with expected i^Pt) 
behavior shown in Fig. 3. 
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q 4 . - - Scenario 2 

. Scenario 3 / 

q g _ Minimum Bias /\ \' 
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P T [GeV] 



(with baryochemical potential which can be in- 
ferred from the baryon transport picture). A dis- 
tinctive characteristic of this picture is that one 
expects qualitatively different vi(j>t) behavior for 
p and 7r and a baryon/meson anomaly limited 
to the intermediate transverse momentum win- 
dow |f-|l]]. Above pt — 2 GeV/c the proton 
and anti-proton v% may keep growing due to flow 
effects, whereas the pion V2, being driven by asym- 
metrically quenched pQCD, may start decreasing 
because strong transverse expansion (if confirmed) 
reduces the energy loss azimuthal anisotropy. The 
comparable contribution of baryons and mesons 
in this region Jl2] , ^9| ^lp^ ] leads to cancelation of 
those effects in the inclusive charged particle ellip- 
tic flow, thus producing a plateau in the 2 < pt < 
5 GeV/c window. At high transverse momenta of 
Pt ~ 6 — 10 GeV/c the pQCD pions dominate the 
charged particle multiplicity, thus leading to a grad- 
ual decrease of V2- The rate at which the high pt 
azimuthal asymmetry vanishes is determined by the 
mean transverse expansion velocity vt, the velocity 
asymmetry (see Fig. 2), and the pQCD fractional 
contribution to the baryon differential pt multiplic- 
ities. The schematic i>2 behavior in this scenario is 
illustrated in Fig. 4. It predicts a detectable differ- 
ence in the moderate to high transverse momentum 
behavior of the elliptic flow for pions and protons. 



FIG. 3. Three possible scenarios of the V2(pt) behavior are 
shown. Error bands are suggestive of the theoretical un- 
certainties within each framework (e.g. the density of the 
medium or the transverse expansion velocity «t). 

1. If experiments confirm an approximate saturation 
of t>2 ~ 0.10 — 0.15 at high pt for both pions and 
protons together with a baryon to meson ratio drop- 
ping back below unity this may suggest that the 
transverse flow might be significantly smaller than 
the values used in the thermal model description 
of relativistic heavy ion collisions. The assumption 
of very large transverse expansion velocities may 
be a caveat since it is not supported by pion inter- 
ferometry as pointed out in Ref. |p9| . In this case 
the energy loss asymmetry plays significant role [^) , 
producing comparable asymmetries for p and n. Al- 
ternatively, appropriate description of particle pro- 
duction could be provided by dissipative relativistic 
hydrodynamics ]30| or inelastic parton cascades |3l]] 
with large partonic cross sections a ~ 10 — 40 mb. 



3. Last, but not least, there are predictions based on 
boosted thermal sources and ideal hydrodynamics. 
The predicted V2(pt) grows continuously as a func- 
tion of the transverse momentum and saturates at 
a maximum value. This was parametrized from hy- 
dro simulations J33| in Ref. || by 



v 2 (p T ) = tanh[p T /(10± 2 GeV/c)] . 



Similarly the baryon/meson ratio (assuming equal 
freeze-out temperatures) is a monotone function of 
Pt and saturates at (p/7r) max = 2 from simple spin 
counting. 



2. Another possibility arises if the pion multiplicity at 
moderate high transverse momenta is dominated 
by the quenched pQCD spectra but the baryon 
production at moderate transverse momenta 2 < 
Pt < 5 GeV/c is controlled by a non-perturbative 
mechanism such as baryon junctions (which ex- 
hibit azimuthal asymmetry and flow) or hydro 



Precision data on identified hadron spectra at high pt 
at RHIC is needed to shed light on which if any of these 
dynamical scenarios predicts the correct flavor depen- 
dence of the azimuthal asymmetry and the proton to pion 
ratio. 
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FIG. 4. Qualitative illustration of the V2(pt) behavior for 
pions, protons and inclusive charged particles within the 
framework of the second scenario is given. A characteris- 
tic difference in the elliptic flow of p and it at moderate high 
Pt is illustrated. Error bands hint at uncertainties related 
to the time for developing a rarefaction wave, flow velocities, 
pQCD vs. hydro vs. string and junction components and 
the baryon/meson flavor composition of charged hadrons. 

IV. CONCLUSIONS 

We have analyzed the effect of possible strong trans- 
verse flow vt ~ 0.6 |l7]|l^,^(| on the transverse azimuthal 
asymmetry generated by jet energy loss. We have shown 
that while the overall magnitude of the quenching is only 
weakly dependent on vt, the azimuthal asymmetry v-i of 
the quenched pQCD component is significantly reduced. 
In order to account for the observed saturation of t>2 of 
charged hadrons, scenario 2 based on (5||nJ is forced into 
predicting an anomalous asymmetry of the baryon com- 
ponent alone. This scenario interprets the anomalous 
baryon to pion ratio in terms of baryon junction dynam- 
ics [|l9] . It differs from scenarios 1 and 3 in predict- 
ing that the pion asymmetry is expected to vanish more 
quickly with increasing px , while the baryon asymmetry 
is expected to be near maximal up to ~ 4 — 7 GeV/c. 
Hydrodynamic and dissipative transport models expect 
a much less dramatic dependence on pt . 
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